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1. Introduction. 

A uniform random fc-dag is an infinite directed graph defined as follows. For each of the integers 
1,2, . .., we pick a random set of fc parents with replacement uniformly from among the smaller non- 
negative integers. This defines an infinite directed acyclic graph (or, dag) with one root (0), and can be 
viewed as a (too) simplistic model of the web, a random recursive circuit (Diaz, Sperna, Spirakis, Toran 
and Tsukiji, 1994, and Tsukiji and Xhafa, 1996), and a generalization of the URRT (uniform random 
recursive tree), which is obtained for fc = 1. All the asymptotic results in the paper remain valid when 
parents are selected without replacement. 

The uniform random fc-dag restricted to vertices 0, 1, . . . , n, is denoted by n or simply U n . 
Indeed, we will take fc = 2 in the main part of the paper, and point out the obvious modifications needed 
when fc > 2 as we proceed. The infinite dag is denoted by Uoo- 

From a given node n, let V n be the collection of paths from node n to the origin. The length of 
path p e V n is L(p). One can consider various path lengths: 

S n = min L(p) , R- = L(P~) , R n = L(P n ) , R+ = L(P+) , L n = max L(p), 

peVn peVn 

where S, R and L are mnemonics for shortest, random, and longest, and P~,P n and P+ are the paths in 
V n , where we follow the parent with the smallest index, the first parent and the parent with the largest 
index, respectively. We can regard R„ and i?+ as greedy approximations of S n and L n respectively. Note 
that, at least in a stochastic sense, 

S n < R n < Rn < Rn < L n . 

The length of the longest path is relevant for the time to compute the value of node n in a random 
recursive circuit, when nodes know their value only when all parents know their value. However, there are 
situations in which node values are determined as soon as one parent or a subset of parents know their 
value — they are called self-time circuits by Codenotti, Gemmell and Simon (1995). For the one-parent 
case, this leads naturally to the study of S n . In networks, in general, shortest paths have been of interest 
almost since they were conceived (Prim, 1957; Dijkstra, 1959). 

It is of interest to study the extreme behavior, as measured by 

max Se , max RJ , max Re , max R~f , max Lp. 

l<i<n l<£<n 1 l<i<n l<£<n 1 l<i<n 

If we replace max by min in these definitions, we obtain the constant 1, and it is therefore more meaningful 
to ask for the exteme minimal behavior as defined by 

min Si , min RJ , min Rf , min Rf , min Lc. 

n/2<i<n n/2<£<n 1 n/2<l<n n/2<l<n n/2<l<n 

So, in all, there are fifteen parameters that could be studied. 
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We take this opportunity to introduce the label process, which will be referred to throughout the 
paper. The label of each parent of n is distributed as \nU\ , with U uniform [0, 1]. An ^-th generation 
ancestor has a label distributed like 

L- • • LLnfiJ U 2 \--- U £ \ G [nUiU 2 ■■■U e -i, nU x U 2 ■ ■ ■ U t ] , 
where the f/j's are i.i.d. uniform [0, 1] random variables. 

The parameter i? n . It is clear that i? n is just the distance from node n in a urrt to its root. In 
particular, R n and its minimal and maximal versions do not depend upon k. We dispense immediately 
with Rn and its extensions because of well-known results on the URRT obtained via the study of branching 
random walks by Devroye (1987) and the equivalence between R n and the number of records in an i.i.d. 
sequence of continuous random variables (see, e.g., Renyi (1962), Pyke (1965), Glick (1975) or Devroye 
(1988)). Only the minimal parameter for i? n requires a gentle intervention. We know that 

1 in probability, 



logn 

for example. Furthermore, 

R n - log n c 



Vlog 



n 



where Af is a standard normal random variable, and — > denotes convergence in distribution. Furthermore, 
an explicit tail bound on R n will be needed further on in the paper. The maximal value of R^ 1 1 < I < n, 
follows immediately from either Devroye (1987) or Pittel (1994). We summarize: 



Theorem 1 . We have 

► 1 in probability, 

logn 

maxKK n R e 

logn 

and 



e in probability, 



lim P \ min Re < 2 i = 1. 

n^oo [ n /2<e<n J 

Finally, for t > log n integer, 

P{R n >t}< cxp (t - logn - tlog(t/logn)) . 



Proof. An outline of proof is needed for the third part and the explicit bound in part four. Let us count 
the number of nodes of index in [l,n/2] that connect directly to the root. This number is 

n/2 

z = ^2^i/e> 
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where £ p is Bernoulli (p). Let A be the event that no node of index in (n/2, n] connects to a node counted 
in Z. This probability is smaller than 

n/2' 



n/2 

<H(l-l/£ + l/(V^i)) 



< exp 




<(L«/2J) 



If the complement of A holds, then clearly, min n / 2 <^< n R(. < 2, and thus, we have shown the third part 
of Theorem 1. Turning to part four, note that R n < min{i : nU\ ■ ■ ■ Ut < 1}, and thus that 

P{Rn >t}< inf EUnUi ■■■U t ) X ) = inf n X (\ + 1)"* = exp (t - logn - tlog(t/ logri)) . □ 
A>0 L J A>0 



Conjecture. For all fifteen parameters, generically denoted by X n , there exist finite constants x 
x(k) > such that 

x in probability. 



logn 



Remark. The limits in the conjecture are denoted by <j, p~, p, p + and A for S n , R^, Rm R-n an d L n , 
respectively. For the minimal and maximal versions of these parameters, we will use the subscripts min 
and max, respectively, as in and cr max , for example. 



Let us briefly survey what is known and provide conjectures in the other cases. 

The parameter L n . Tsukiji and Xhafa (1996) showed that A max = ke. The Chernoff large deviation 
bound shows that A is at most the largest solution x of 



x 



e-' = l, (1) 



and thus A < A max . We believe that A is indeed given by (1) based on arguments not unlike the proof of 
Theorem 2 below. We have no guess at this point about the value of A m ; n . 



The parameter i?+ . In the label process, the parent's index is approximately distributed as n max(C7i , ... ,11^), 
where the C/j's are i.i.d. uniform [0, 1] random variables. If U, as elsewhere in this paper, is uniform [0, 1], 
then the parent's index is thus roughly like nU 1 ^. By renewal theory, this implies that 

> k A = p + in probability. 

log n 
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Chernoff's large deviation bound show that Pmax ls & t most the unique solution x of (2) that is above k: 

e l - k = 1. (2) 



ke\ 



x 



We believe that the solution of (2) yields p^ax • Applying Chcrnoff to the other tail shows that p^in ^ s 
at least the other solution of (2), as (2) has two solutions, one below k and one above k. Furthermore, 
we believe that the two solutions of (2) yield p+ in and p+ ax . 

For k = 2, the parameter R^ is intimately linked to the random binary search tree, which can 
be grown incrementally by a well-known process described as follows: given an n-nodc random binary 
search tree, sample one of its n + 1 external nodes uniformly at random, replace it by node n + 1, and 
continue. The parent of that node is either its neighbor (in the total ordering) to the left or its neighbor 
to the right, and in fact, it is the neighbor added last to the tree. But the labels (times of insertion) of 
the neighbors arc uniformly drawn without replacement from {1, . . . , n}, and are thus roughly distributed 
as nU, so that the parent of n + 1 is roughly distributed as nVU, because the maximum of two i.i.d. 
uniform [0, 1] random variables is distributed as \f\J . With this in mind, max 1 <^< n i?^" is the height of 
the random binary search tree, i?+ is the depth (distance to the root) of the node of label n (the n-th 
node inserted), and min n / 2 <^< ri Rg is very roughly the shortest distance from leaf to root, or fill-up level. 
These quantities behave in probability as described above, as shown by Devroye (1986, 1987), and this 
explains the values p+ ax = 4.31107 . . ., p+ = 2 and p+ in = 0.3733 . . .. 



The parameter R~ . Arguing as above, the parent's index is approximately distributed as n mm(Lq , . . . , U k ). 
By a property of the uniform (or exponential) distribution, using a sequence of i.i.d. exponential random 
variables E\ , E<i , . . . , we have this distributional identity: 

nmin(£7i, . . . , U k ) = nUiU^ 2 • • • U^ k = exp I logn — -4- ] ■ 

V j=i 3 / 

Renewal theory easily gives the law of large numbers and central limit theorem for R~. For example, 

■ — i- = f p~ in probability, 



logn H k 

where — X}j=l(Vi) i s the fc-th harmonic number. Using large deviation bounds similar to the ones 
used below in showing part of Theorem 2, one gets that 



lim P < max R . > (x + e) log n > = 

n^oo [l<e< n 1 J 

for all e > 0, where x is the solution greater than l/H^ of 

k 

l + f(x)=x^og(l + f(x)/j), 

3=1 

and f(x) > is implicitly defined by 

k 1 1 

E- = -, x > l/H k . 
■^3 + f{?) * 
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These equations follow from the obvious Chernoff bound. We conjecture that p max equals this upper 
bound, but a rigorous proof that p max is indeed as described above is not given in this paper. 



The parameter S n . The most important parameter for computer scientists and combinatorialists is the 
one in which graph distances are defined by shortest paths, and this leads to the study of S n . That was 
the original motivation of the paper, and we will settle first order asymptotics in this paper. Theorem 1 
implies, for example, that with probability tending to one, 

min Sg < 2. 

n/2<£<n 

So we turn to a and er max : 



Theorem 2. Assume k>2. Then a = <7 max , where a is given by the solution x e (0, 1) of 

(£) V = l. (3) 
[Note that ip is indeed an increasing function on (0, 1) J 

Observe that Theorem 2 does not extend to k = 1, because in that case, S n = R n = L n , and 
similarly for the maximal versions of these parameters, in view of the equivalence with the urrt. Thus, 
S n / log n — ► 1 and max 1 <^< ra Sg/ logn — > e in probability. 

The following is a table of constants in the Conjecture for k = 2. The constants involving a 
(top row) are obtained in this paper, while those involving p (third row) are covered by Theorem 1 . The 
constants p~ and p + follow from ordinary renewal theory. The zeroes in the table follow from Theorem 
1. Finally, A m ax is due to Tsukiji and Xhafa (1996). There are thus four conjectured constants, which 
happen to be one-sided bounds (p max , P mm j Pmax ■> an d one unknown constant, A mm . 



^min ® "max 
/'min P /'max 
Pmin P Pmax 

PmmP + Pt 



min f f max 
X max 



A m in A A 



0.3733. . . 0.3733. . 

0.6666... (= 2/3) 1.6737.. 

1 2.7182.. 

0.3733... 2 4.3110.. 

? 4.3110... 5.4365.. 



(= e) 
(= 2e) 



The following is a table of a, p and p max for different numbers of parents k. 
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k 


a 


P~ 


Pmax 


2 


0.3733... 


0.6666. . . 


1.6737... 


3 


0.3040. . . 


0.5454. . . 


1.3025... 


4 


0.2708. . . 


0.48 


1.1060... 


5 


0.2503... 


0.4379. . . 


0.9818... 


6 


0.2361... 


0.4081... 


0.8951... 


7 


0.2254. . . 


0.3856... 


0.8305. . . 


8 


0.2170... 


0.3679. . . 


0.7800. . . 


9 


0.2102... 


0.3534. . . 


0.7393. . . 


10 


0.2045. . . 


0.3414... 


0.7057. . . 


11 


0.1996... 


0.3311... 


0.6773. . . 


12 


0.1954... 


0.3222... 


0.6531... 


13 


0.1916... 


0.3144. . . 


0.6318... 


14 


0.1883... 


0.3075. . . 


0.6132... 


15 


0.1854... 


0.3013... 


0.5966. . . 


16 


0.1827... 


0.2957. . . 


0.5816... 


17 


0.1802... 


0.2907. . . 


0.5683. . . 


18 


0.1780... 


0.2861... 


0.5560. . . 


19 


0.1760... 


0.2818... 


0.5448. . . 


20 


0.1740... 


0.2779. . . 


0.5346. . . 


21 


0.1723... 


0.2743. . . 


0.5251... 


22 


0.1706... 


0.2709. . . 


0.5164... 


23 


0.1691... 


0.2677. . . 


0.5083. . . 


24 


0.1676... 


0.2648. . . 


0.5007. . . 


25 


0.1663... 


0.2620. . . 


0.4936. . . 


26 


0.1650... 


0.2594. . . 


0.4868. . . 


27 


0.1638... 


0.2569. . . 


0.4805. . . 


28 


0.1626... 


0.2546. . . 


0.4747. . . 


29 


0.1615... 


0.2524. . . 


0.4690. . . 


30 


0.1604... 


0.2503. . . 


0.4638. . . 


35 


0.1559... 


0.2411... 


0.4409. . . 


40 


0.1521... 


0.2337. . . 


0.4225. . . 


45 


0.1490... 


0.2275. . . 


0.4074. . . 


50 


0.1463... 


0.2222. . . 


0.3946. . . 
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2. The shortest path length S 1 



We will establish Theorem 2 in two parts. First we show that for all e > 0, 

lim P{S n < (1 - e)crlogn} = 0, (4) 

and then that 

lim P I max St > (1 + eWlogn I = 0. (5) 
We only consider the case fc = 2 since the case k > 2 follows quite easily. 



Lemma 1. Let G a be gamma(a) , with a > 1. Then 



X r(a) 

and 



< ^—^T,x>a-1, 



P{G a < x} 1 



Proof. The gamma density is f{y) = y a e y /T(a). It is log-concave for a > 1, and thus, a first-term 
Taylor series bound yields the inequality 

f(y) < f(x)e ( y- x ^ ogf y^ = f{x)e^- x ^ {a -^l x - l \ 

Integrating the upper bound out over [x, oo) or (— oo,x] then immediately yields the results. □ 

From node n, we can consider the index of the first of the 2 £ i-th level ancestors, which is 
distributed as 

L- • • [[nUi\ U 2 \ ■ ■ ■ U e \ > nUiU 2 ■■■U l -(. = nexp(-G £ ) - t, 

£ 

where = denotes equality in distribution, and is gamma(l). If these indices are I\, . . . ,I 2 e, then we 
have 

P{S n < 1} = P < min h = \ 
[l<i<2^ I 

< 2^P{7i = 0} 

< 2 e P{ncxv(-G e ) -I < 0} 
= 2 f P{G e > \og{n/ 1)} 

„ <? 3 / 2 (21og(n))V lo sW 



Set I = [t\ogn\ for t e (0, 1), and note that the upper bound is 

e (io g 3 /2 (n) ) x {<p{t))^ n , 

where (p(t) = (2e/t) t /e is as in (3). We have if (a) = 1 for a = 0.3733 . . .. Thus, we have shown (4): for 
all e > 0, 

P{Sn < (a-e)logn} = o(l). 

Although we will not need it directly, we will also deal with the upper bound on S n . This can 
be done in a number of ways, but the shortest route is perhaps via the great-grandparent strategy that 
jumps I generations at a time, where I now is a large but fixed integer. We denote this by i-GGP. We 
associate with each node n two independent uniform [0, 1] integers U and V and let the parent labels be 
\nU\ and [nV\ . Let A n be the event that any of the 2^ ancestors of node n conicidc. It is clear that 
P{^4n} — *• as n — > oo. As an ancestor label is described by 



we define 



L- • • llnU^Uil ■ ■ ■ U e \ < n\J x U 2 ■■■U l ^ nexp(-G £ ), 



min TT U e 



where V is the collection of all paths of length I above node n, and each pe? consists of edges e that 
each have an independent uniform random variable associated with it. If e > and n is greater than 
some n e , then the £-GGP gives with probability greater than 1 — e a node with label less than Zpa. Define 

^ (e) =min(^,6) 

where b is chosen such that V{Zi > b} = e. As long as the label stays above n e , one can dominate the 
labels in the £-GGP by multiplying n with successive independent copies of Let T n be the number 

of steps until the label in ^-GGP reaches n e or less. Renewal theory shows that with probability tending 
to one, 

(1 + e) logn 

"- E {-.o g ( z <«')}' 

because the £-GGP takes I steps at a time, and because a node with label n e is not further than n e away 
from the origin, we see that with probability tending to one, 

1(1 + e) logn 



S n < n e + 



< 



E{-log«)} 
+ 2e) logn 



Uniform integrability implies that 

lim E | - log (Z^) } = E {- log (Z t )} . 

Therefore, for any (new, fresh) e > and £ > 1, with probability going to one, 

+ e) logn 



E{-log(^)}' 
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Observe that 

- log ( z e) c 1 v-^ j? 
1 — = 1 ma ^ 1^ Eu > 

eEp 

where the E u are i.i.d. exponential random variables. From the theory of branching random walks, it is 
easy to verify (see, e.g., Biggins (1977), or Devroye (1986, 1987)) that, as I — > oo, 



1 1 
- max > E u -» - 

in probability. Thus, 



p&V , 
eEp 



liminf -E{log(^) }> l ; 
£->oo £ CT 

and thus, by choosing £ large enough, we see that with probability tending to one, 

S n < (l + 2e)crlogn. 

This concludes the proof of the first part of Theorem 2. 

The next section requires an explicit rate of convergence. To this end, still restricting ourselves 
to k — 2 only, let zj^, 4^' ... be i.i.d. copies of zjj € \ and note that, 

T n < min jt : nZf^ ■ ■ ■ Z^J < 1 j 

= min {t : log (l/4j) + ' ' ' + lo S (V^,?) > lo S n } ■ 

Set n = E | log (l/^) }■ Thcn > assuming 5* e (0, 1/2) and <5 e (<5*, 2(5*) such that m = (l//z + S) logn 
is integer- valued, 

P{T„ > m} < P {log (l/4i) + ' ' ' + log (l/4t) < lo & n } 

= P {log (l/4,i) + • • • + log (l/4i) — m M < -5/xlogn} . 

Let p > 2 be a fixed number. Rosenthal's inequality (Rosenthal, 1970, Fuk and Nagaev, 1971, 
see also Petrov, 1975) states that there is a constant C p with the following property. If {X n ,n > 1} is a 

sequence of centered and independent random variables, and if Y n = X\ H h X n , and if E{| X n | p } < oo 

for all n, then 

H\Yn\ P } < C p ^ H\Xj\ P } + (V{Yn}) P/2 

For i.i.d. random variables with X\ = X, we have 



E{|Fn| p } < C p (nE{|X|P} + n P/ 2 (E{X 2 }) j < 2C p max(n,rW 2 )E{|AT}. 
Applied to our situation with p = 4, using Markov's inequality, we have 

P{T n > m} < (<5 M logn)- 4 E |(log (l/4,i) + ' ' ' + lo S (v4,l) ~ TO ) 

< 2C 4 (^logn)- 4 m 2 E | log (l/4 £) ) - M 

< C(logn)- 2 ^" 4 , 
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where C depends upon e and I only. The remainder of the argument involving an appropriate choice of 
I remains valid, and we can conclude that for any e > 0, 

P{S n > (<r + e) logn} = O (l/ log 2 n) , (6) 

with room to spare. 



3. The maximal shortest path length 



The purpose of this section is to show (5). We let a be as in the first part of the proof, and let 
e > be arbitrary. Fix n large enough. From (6), 

E {\{j : n/2 < j < n, Sj > (a + e) logn}| } = O 



log 2 n 



and thus P{A(n)} = O (^^) , where 



A(n) 



dcf 



{j : n/2 < j < n, Sj > (cr + e) \ogn}\ > J}] ■ 



If we take an incremental view of the process of adding edges, then a node with index in [n, 2n] selects a 
parent of depth < (cr + e) logn and index > n/2 with probability > 1/8 if A(n) fails to hold. It is this 
observation that will allow us to uniformly bound all depths by something close to (cr + e) logn. 

Consider the indices in dyadic groups, {2 r - 1 + 1, . . . , 2 r }, r > 1. We recall from a comparison 
with the urrt, that S n < R n and thus that maxi<j< n Sj < maxi<j<„ Rj, and that (see Theorem 1) 



P j max R j > 2elognl < n~ 2elog ( 2 ) 
{l<j<n J J " 



< n 



Thus, for 7 > small enough, 

P 

It remains to show that 
Consider the event 



i<j<L" 7 J 



} 



max Sj > (cr + e) logn \ = 0{n 37 ) = o(l). 



P < max Sj > (cr + e) log n 



n1<i<n 



(!)• 



B(r)= |J A{2 S ), 



r <s<r 



where r' is the largest integer such that 2 r < n 7 . Clearly, P{B(r)} = 0(l/r') = 0(1/ logn). On the 
complement, (B(r)) c , intersected with max 1 <j<^ ri 7j Sj < (cr + e) log(n) , we look at the process started 
at a node m < n and assume that its index m is in |2 r + 1, . . . , 2 r+1 }. That process is looked at as a 
binary tree of consecutive parents, and will be cut off at height h = [10 log log nj . There may be duplicate 
parents (in which case the tree degenerates to a dag) , so we need to be a bit careful. If any parent in the 
tree is selected with index < 2 r < n 7 , then S m < (cr + e) logn + h, and thus, we can assume that in this 
"tree" any node j selects its parent uniformly in the range (2 r , j). At any stage, by our assumption, the 
probability of picking a parent i having Si < (a + e) logn is at least 1/8 (and this is why we needed the 
dyadic trick, so that we can make this statement regardless of the choice of i within the range (2 r , n]). 
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We claim that this "tree" has at least 2 fe ~ 1 leaves or reaches [1, 2 r ] with overwhelming probability. To 
see this, note that a node j in it picks a node already selected with probability not exceeding 2 h /j. But 
the index j is stochastically larger than 

X h d ^ [■■■[[mU 1 \U 2 \---U h \ 

by our remarks about the labeling process. The probability that there are in fact at least two such 
unwanted parent selections (but none of them less than n 7 ) in that "tree" is not more than 

^ >< ^ * ^ >< ^ {j^IW]} (?) 

We have 

■dt 



/■OO 

= / P{X h <l/t}dt 

Jo 



10 

< / ¥{mUi ■■■U h <h + 1/t} dt 
Jo 



-I 

Jo 

= / / " r(h\ d y dt 

o° y h-l e -y f _^ | 



777 

p ^x+i7t <Gfc} * 

-.l/n 7 /-oo yh-l e -y 



T(h) " V (me-v-h) 

^ / fk^ dy + n ^ 

Jo T{h)m J k 



dy 

,-v 

dy 



\og(m/2h) r(/l) 

2(log(n)) fe n-^logH)^ 1 ^ 

< — 1 — (for n large enough, by lemma 1) 

mh\ T{h)m 



O (n°W/m) = O (m 



-l+o(l) 



Thus, our probability (7) is not larger than O If there is only one unwanted parent selection 

and we avoid indices below n 7 , and considering that the first parent selection at the root node is always 
good, we see that at least half of the 2 h potential leaves are in fact realized. Each of these leaves makes two 
independent parent selections. The probability that all these leaves avoid parents j with Sj < (c + e) logn 
is at most (7/8) 2 = o(n~ 2 ). If there is a connection, however, to such a parent of low depth, then 
the root has shortest path length at most h+1 more than (a + e)logn. Hence, if £ m is the event 
[S m > (o- + e) log" + h + 1], then 



p\£ m n((B(r)) c n 



max Sn < (a + e) log n 

i<i<L« 7 J . 



|=0 (m" 2 ^ 1 )) . 
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Thus 



P<J max Si > (er + e)logra + /i+ll 
\n't<j<n J J 



— P {u^> n7 £ m } 

<p( max Sj > (cr + e)lognl + P{B(r)} + V m _2+ °W 

U<i<L«TJ J m >„ 7 



= 0(n" 37 ) + 0(l/r') + n-^ 1 ) 
= 0(1/ log n). 

This concludes the proof of the theorem. □ 



4. Bibliographic remarks and possible extensions. 



The study of the URRT goes back as far as Na and Rapoport (1970) and Meir and Moon (1978). 
Single nonuniform parent selections have been considered as early as 1987 by Szymahski. Szymanski 
(1987) showed that if a parent is selected with probability proportional to its degree, then with high 
probability there is a node of degree ^l(^/n). This is nothing but the preferential attachment model of 
Barabasi and Albert (see Albert, Barabasi and Jeong, 1999, or Albert and Barabasi, 1999), which for a 
single parent is a special case of the linear recursive trees or port (plane-oriented recursive tree). For 
this model, the parameter R n was studied by Mahmoud (1992a), and the height by Pittel (1994) and 
Biggins and Grey (1997), and in a rather general setting by Broutin and Devroye (2006): the height is 
in probability (1.7956 . . . + o(l)) logn. The profile (number of nodes at each depth level) was studied by 
Hwang (2005, 2007) and Sulzbach (2008). 



One can ask the questions studied in the present paper for these more general models. 



Various aspects of URRT's besides the depth and height have been studied by many researchers. 
These include the degrees of the nodes, the profile, sizes of certain subtrees of certain nodes, the number 
of leaves, and so forth. Surveys and references can be found in the book by Mahmoud (1992b) or the 
paper by Devroye (1998). Specific early papers include Timofeev (1984), Gastwirth (1997), Dondajewski 
and Szymanski (1982), Mahmoud (1991), Mahmoud and Smythe (1991), Smythe and Mahmoud (1994), 
Szymahski (1990), and the most recent contributions include Fuchs, Hwang and Neiningcr (2006), and 
Drmota, Janson and Neininger (2008). One may wonder how the profiles behave for uniform random 
fc-dags. 
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